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In the framework of the classical effective Lorentz invariance violating (LIV) model of Myers- 
Pospelov, we present a complete calculation of the synchrotron radiation produced by a circularly 
moving charge in the rest frame of the model. Within the full far-field approximation we com- 
pute exact expressions for the electric and magnetic fields, the angular distribution of the power 
spectrum and the total emitted power in the m-th harmonic. We also perform an expansion of 
the latter quantity in terms of the electromagnetic LIV parameter and calculate the average degree 
.^j- ■ of circular polarization to first order in such a parameter. In both cases we find, under adequate 

' circumstances, the appearance of rather unexpected and large amplifying factors, which go together 

, with the otherwise negligible naive expansion parameter. This opens up the possibility of selecting 

■ astrophysical sources where these amplifying factors are important, to explore further constraints 

^ ' imposed upon the LIV parameters by synchrotron radiation measurements. 

CD 



PACS numbers: 11.30.Cp, 41.60.Ap, 03.50.Kk, 95.30.Gv 



I. INTRODUCTION 

£^ ' The possibility of exploring Lorentz invariance violations (LIV) at the highest possible energies has been further 
, motivated with the proposal of Ref. Q that photons may acquire modified energy-dependent dispersion relations 
t-H ' (MDR). The effects on the photon velocity could be substantially amplified when looking both at cosmological and 
high energy sources (gamma ray burts, for example) to allow for an observational detection of such minute corrections. 
This idea has been further generalized to massive spin 1/2 particles, in such a way that the typical form of such modified 
dispersion relations is 

O ; 

P^; c 2 (fc) = fc 2 ±^, E 2 (p)=p 2 + m 2 + VR , L ^, (1) 

Oh! respectively. Here M is a scale which signals the onset of LIV and is usually associated with the Planck mass Mp. 
This suggests that quantum gravity could be responsible for such a breaking of Lorentz covariance. This proposal 
has been further explored in Loop Quantum Gravity, where constructions of increasing degree of sophistication have 
managed to produce modified effective actions for photons and fcrmions encoding the relations |QJ. Basically, the 
starting point of these approaches has been the well defined Hamiltonian operators of the quantum theory together 

5^ 1 with an heuristical characterization of the semiclassicalground state 0, 0, 0, HJ . String theory has also provided a 
possible connection between quantum gravity and LIV |al3- 

The MDR JTJ generically arise from Lorentz-violating dimension five operators in the Lagrangian, upon which 
we focus in the following. On the phenomenological side the main emphasis has been set on the determination of 
observational bounds for the LIV parameters appearing in the MDR. A partial list of references is given in 0, [jj . 
Additionally, effective field theories, not necessarily related to quantum gravity, have been put forward to analyze 
the possible consequences of LIV [HJ, Ell- In this work we restrict ourselves to classical theories, and therefore we 
do not address the fine-tuning problems arising from quantum corrections associated with the coexistence of space 
granularity and preferred reference frames [T^ . 

Here we are concerned with the description of synchrotron radiation coming from cosmological sources, which 
have already been used to set bounds upon the parameters , VR m O @- Such bounds are based on a set of 
very reasonable assumptions on how some of the standard results of synchrotron radiation extend to the Lorentz 
non-invariant situation. This certainly implies some dynamical assumptions, besides the purely kinematical ones 
embodied in More specifically, we examine these assumptions in the light of a particular model, which we 

choose to be the classical version of the Myers-Pospelov effective theory (MPET), which parameterizes LIV using 
dimension five operators Besides, a complete calculation of synchrotron radiation in the context of this model 
is presented. This constitutes an interesting problem on its own whose resolution will subsequently allow the use of 
additional observational information to put bounds upon the correction parameters. We have in mind, for example, the 
polarization measurements from cosmological sources. The case of gamma ray bursts has recently become increasingly 
relevant 0, although it is still at a controversial stage Q. It would also be interesting to compare the results from the 
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MPET with those arising from the Gambini-PullinQ, as well as the Ellis et al.Q electrodynamics. We postpone these 
questions, together with most of the details in the MPET calculation, for separate publications [1 3 . Our calculation 
rest heavily upon the work by Schwinger et al. on synchrotron radiation, reported in Refs. [Tfil. ll7l . ll8l ]. A partial list 
of previous studies of electrodynamics incorporating LIV via dimension three and four operators is given in Ref. |l9j| . 



II. MYERS-POSPELOV ELECTRODYNAMICS 



The Myers-Pospelov approach is based on an effective field theory that describes Lorentz violations generated by 
dimension five operators parameterized by the velocity of a preferred reference frame, not taken as a dynamical 
field. These operators are assumed to be suppressed by a large inverse mass factor, which we take to be Mp 1 , and 
hence can be considered as perturbations at the classical level. In this framework we analyze the electromagnetic 
radiation of a classical charged particle. To avoid unnecessary complications we restrict the discussion to the case 
V 1 - 1 = (1, 0), enough to show the main consequences of the Lorentz violations. We take the unperturbed velocity of 
light c in vacuum equal to one. 

The dynamics of a classical charged particle can be obtained from the geometrical optics limit of a scalar charged 
field. In this case the Myers-Pospelov action is 



5' 



MP 



dpip'dfip - n*tp*tp + irj<p* (V ■ dy if 



(2) 



with V ■ d = d^. In momentum space (ip(x) = (fio expi(p°t — p • x)) and in the reference frame where V a — (1,0), 
the dispersion relation becomes 



(p°) 2 + v(p°) 



(3) 



with fj — —r)/Mp, where rj is the dimensionless constant employed in the parameterization of Jacobson et al. 9]. The 
Eq. (|3J) is an exact relation in fj. From here we obtain the Hamiltonian for a massive particle, to second order in fj 



2\l/2 



2\3/2 



£3 \ 



(4) 



In what follows we consider the interaction of a particle of mass /i and charge q with a static magnetic field 
{4> = 0, A(xi)) via the standard minimal coupling. From the corresponding Hamiltonian equations we obtain the 
acceleration 



i- = |( l -■-,-,£ + ^) { y <B). 



(5) 



with E = H(p — qA) being the energy of the particle. As in the usual case, this means that: (i) the magnitude |v| is 
constant and (ii) the projection of the orbit in a plane orthogonal to B = V x A is circular, with a Larmor frequency 



LOq 



E 



1 



-fjE- 
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fj 2 E 2 



(6) 



Here we consider only the circular motion in a plane orthogonal to B and we take the z-axis in the direction of the 
constant magnetic field. Using the standard definition (3 = |v|, the solution to the equations of motion can be written 
as 



8 

r(i)= — (cos wot, sincjot, 0) . 



(7) 



which identifies the Larmor radius of the orbit as R = (3/ujq. The modified relation defining in terms of the particle 
energy is 



l-(3 2 



El 

E 2 



1 



fjE 3 15 f] 2 E A 



Let us remark that we maintain the standard definition of the Lorentz factor as 7 
by the corresponding charged particle is 

j(t,r) =^ 3 (r-r(t))r(i), 



!3 2 ) 1/2 . 



(8) 

The current produced 
(9) 
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where r(t) is given by Eq.JTJ). 

The Myers-Pospelov action for the electromagnetic field is 



Smp = J d 4 x -\f^ v - 4nj»A„ + £(V a F aS ) (V ■ d){VpF^ s ) 



(10) 



where £ = £/Mp, with £ being a dimensionlcss parameter. As usual we have — d^A v — dyA^, j M = (p,j) and 
we are using the conventions in Jackson |2dj| . Clearly this is a gauge theory and implies the conservation of the 
electromagnetic current j M . The definition of F^ leads to the standard homogeneous Maxwell equations. The source 
dependent equations of motion obtained from the action Hl()(l . written in the rest frame V a — (1,0), are 



<9E 



d 



V • E = 47rp, — — + V x B + £— -V x E + — - = 4?rj 



dt 



dt 



dB 



dt 



(11) 



Introducing the standard vector potential field A^ = (cf>, A) and using the radiation gauge V • A = we are left with 

8 2 A ~ 8 2 A ( 1 \ 

V 2 A + 2£Vx — =4tt j-V— V-j =47rj T . (12) 



dt 2 * dt 2 

Let us recall that in this gauge the electric and magnetic fields in the radiation approximation reduce to 

BA 

E = -^, B = VxA. (13) 
The energy momentum tensor T^ v for this modified electrodynamics is given by 

(14) 



S=— E x B — E x 

4tt 4tt dt 



which satisfies the usual conservation equation dtT 00 + V • S = outside the sources. These expressions are exact in 
I 

To solve the equation of motion for A it is convenient to work in momentum space, where Eq. I|12|) becomes 



(-Lu 2 + k 2 -2ifw 2 kx) A(w,k) =47rj T (w,k), 



(15) 



with jT(w,k) = j — k (^k • jj and k = k/fc with k = |k|. The notation we employ for the different types of Fourier 

transforms is via the corresponding arguments. For example, if F(t,r) denotes the function in space-time, F(u),r) 
denotes the Fourier transformed function to frequency space, while F(u>,\s.) denotes the Fourier transformed function 
to frequency and momentum spaces. Equation l|15|) can be diagonalized using a circular polarization basis yielding 
radiation fields of definite helicity 



A>,k) 



l r 



(t.A) 



k + iXk x A 



with A = ±1, leading to 



[-uj 2 + k 2 - A2£> 2 fc) A A (w, k) = 47rj£(w, k), 



where 



1 



1 



j£(w,k) = - j T + zAkx j T = - j- k-j k + Ukxj = j A (w,k). 



The Eqs. i|17|) imply the following dispersion relations for the propagating fields 

u\ = — ^. 

1 + 2\£k 

To first order in £, they have the same form as those corresponding to the photon in Eq. (JJJ. 



(16) 



(17) 



(18) 



(19) 
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Equation (|17fl can be solved in terms of the retarded Green functions with definite helicity 

f 



G? ei H k ) = 
It is convenient to calculate 



k 2 - X2& 2 k 



(20) 

uj — >LU-\-ie 



G A et H r-r') = / -^ e * (r-r') G A ( k) = 1 ^f) e in(A^ (21) 
re * 7 (2tt) 3 47ri?vTTi2 



rf 3 k 

where i? = |r — r'|. Here we introduce the polarization-dependent refraction index n(Xz) 



n(Xz) = + z 2 + Xz, z = £uj. (22) 

In this way, the fields A A in Eq. (|17|l have well defined phase velocities v\ — l/n(Xz) and the situation can be 
described as the propagation of photons in a dispersive birefringent medium. 

The Green functions H21|l determine the corresponding potentials in the far field approximation with the standard 
replacements l/R ~ l/\r\ = 1/r in the denominator and R ~ r — n ■ r' in the phase, with n = r/r been the direction 
of observation. Using l|21|l we obtain 

A>,r) = 47r J <Ft> GUu,T-T>)&(uJ)=±-10L e *»M<>r J d Ve-«^W y j^(,,r'), (23) 

and thus we finally get 

A>,r) = e MA*Vr-A ( kA) (24) 

in the radiation approximation. Let us emphasize that the momenta k A = n(Az)cjfi are fixed in terms of the frequency 
and the direction of observation. As usual, the integration over d 3 r' has been conveniently written as the Fourier 
transform in momentum space jj>(w, k A ) of the function (u>, r') in coordinate space. In the last term of (|24|) we have 
used the relation l|18f) . The full vector potential is given by the superposition 

A(w,r) = - J—y V n(A*)e*™^'j A Hk A ). (25) 
The magnetic and electric fields satisfy 



B(w,r) = VI + z 2 n x E(w,r) - izE(w,r). (26) 

At this point it is interesting to observe that the above relation, together with the expression for the Poynting vector 
in H14|). yields the z-exact, manifestly positive-definite result 

S = ^v/TT^ E» • EH = ^ + E*H • EH. (27) 

This reproduces the standard expression for the Poynting vector in a medium with refraction index n by the substi- 
tution (n(z) + n(—z))/2 — > n. 

The use of Eas. i|13|l and i|25[) leads to the following electric and magnetic fields in the radiation approximation 

BHr) = ~^L= T Xn 2 (Xz)e m( - Xz ^ r j x (uj, k A ), (28) 

E(w,r) = --=^ £ «(Az) e m ^^j A ( w ,k A ). (29) 

Finally, using for example the methods in Ref. we arrive at the following general expression for the angular 

distribution of the power spectrum in the MPET 



d 2 P(T) 1 



/DC 
dr Yl n 2 (Xz)f k (T + T/2,k x ) P, A > (T - r/2, k A ) (30) 

-CO \ i -i 



dwdn \z 2 ,T~^J , v n 

The projector P£ r can be read from the relation among the components of the polarized and total currents given in 
Eq. Q1H[). The total current is written here in the representation j r (t, k). 
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III. SYNCHROTRON RADIATION 



We now consider the radiation of a charged particle moving in a circular orbit, in a plane perpendicular to a 
constant magnetic field B parallel to the z-axis, according to Eq.Q. The direction of observation is given by 
n = r/r = (sin 6, 0,cos(9) and the products of currents that appear in Ea. (|3U[l lead us to consider the following 
combinations of the velocity v(t)=r(t), according to Expression (JSJ, 

v(T + r/2) • v(T - r/2) = q 2 f3 2 cos uj q t, (31) 
n ■ v(T + r/2) x v(T - r/2) = -g 2 /3 2 cos sin w r. (32) 

The angular distribution of the radiated power spectrum is given by the average over the macroscopic time T of 
Ea. (|30J) . As in the usual case, it is convenient to write the relevant quantities as a sum over the contribution of the 
harmonics m = 0, 1, 2, By doing this we have 

12 



d 2 P(T) \ ^ sr, ( , dP m , x 

' X=±l m=0 



with cj to = mujQ, z rn = £w m and 



— ^ = jp . [Am/3n(Az m )J^(W Aro ) + m cot 9 J m (W Xm )f . (34) 

Here ( ) denotes the average over T, and J m , are the Bessel functions of order m and their derivatives respectively. 
The argument of the Bessel functions is 

W\ m = m n(\z m )(3 sin 0. (35) 

Let us remark that the unpolarized angular distribution of the power obtained from Eq. (|34|l coincides with the 
standard result given in Eq. (38.42) of Ref.pt in the limit £ = 0. From ^ we obtain that the angular distribution 
of the averaged power is anisotropic to first order in £. We also calculate the total averaged and integrated power 
radiated into the m-th harmonic 

q 2 mujQ \ - 



P m = ^ }2 n (^ z m) [2/3V 2 ' m (2m n(\z m )l3) 

IK y/l + Z m A = ±1 

/•2mn(Xz m )/3 

- [n- 2 (Xz m ) - (3 2 ] / dx J 2m (x) 







(36) 



which clearly indicates the contribution of each polarization. The above result is exact in z m and the parity-violating 
contribution vanishes after the angular integration. 

Additional information regarding the LIV parameter £ can be obtained from polarization measurements of the 
incoming radiation. To this end we have calculated the exact expression for the m-th contribution to the electric field 
in the radiation approximation = E;^(w m , k A ) 

= -Ay ; r *P"rnn(\z m )„ in(Xzm)Umr / dJ m {W Xm ) + x J m (W Xm ) \ 
r yfl + z^ V dWxm tan6< J 

where e A denotes the circularly polarized basis 



1 

V2" 



e± = —f= (e|| ± ie±) , ey = (0,1,0), ex = (-cos 61,0, sin 0), (38) 



in the chosen reference frame. The expression l|37|) allows for the calculation of the corresponding Stokes parameters. 

There are two ways in which the photon LIV parameter £ enters the modified expressions for the physical quantities: 
(i) one is via over-all multiplicative factors such as 1/Vl + z 2 or n(Xz), where the dependence upon A£w makes these 
contributions negligible for the range of observed frequencies, (ii) the second possibility is through the dependence 
upon the variable Wx m in the arguments of the Bessel functions. We have found that this case provides some 
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TABLE I: Typical parameters of some relevant astrophysical objects 
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possibilities of producing observable effects because of the appearance of additional amplifying factors, as we show 
below. Now, the contributions to case (ii) are governed by the far field expansion of the Green function phase 

n{Xz)uj |r - r'| ~ ur ^1 - + )£u - A^^- + ^ ( i- ) ), (39) 

where the term proportional to (r'/r) 2 is always neglected in the radiation approximation. Consistency demands that 
one had better make sure that terms proportional to the LIV parameter £ are larger that the neglected one in order to 
properly include them in the argument of the Bessel functions. To this end, in Table^we make a rough estimation of 
the relevant parameters corresponding to some observed cosmological objects. The notation is: r [l.y.] is the distance 
of the object to the earth, 7 is the Lorentz factor of the charged particles, B [Gauss] is the average magnetic field 
producing the synchrotron radiation, to [GeV] is the maximum observed frequency, ujq [GeV] is the Larmor frequency 
and r' = R is an estimation of the size of the radiating source, where R is the Larmor radius. Our general results are 
presented in the full far-field approximation. 

As indicated in Tabled the radiation of interest is dominated by very high harmonics m » 1. Th us, it is convenient 
to use the large m approximation for functions of the type J m (mx) together with their derivatives 16] . In the regime 
1 - n 2 f3 2 > we find 



q 2 mc 1 \( 3 \ 2/3 , 

Pxm = 7n n T~ — JT\ — T 1 ^ — / dxK 5/3 (x) 




(40) 

-'"-Ac/ \'"Ac/ I 

which introduces the cut-off harmonic 

m Xc = \ [l-(3 2 n 2 (Xz m )y 3/ \ (41) 

such that for m > fh\ c the power decreases as P\ m ~ e -™/™Ac This is basically the result of Ref. 01 extended to po- 
larized radiation. Here K p / q denote the Bessel functions of fractional index. Within the same large-m approximation, 
the integrated power in the m-th harmonic can be expanded to second order in £ yielding 

^ ^ K (-)-- 2 K v J™)+2(^)\ v J^)}, ,42, 



VSttR-j 2 I m \m c J j 2 \m c J \ 7 



in,. 



where m c — 37 3 /2 and n(x) — x J°° dy K 5 / 3 (y) is the bremstrahlung function [2l|. Let us emphasize the appearance 
of the combination £uim/j = ^(uj/Mp)(m/j) as the expansion parameter in l|42l) . From Table |U we see that this 
combination is not necessarily a small number, which signals the possibility that such corrections might be relevant 
in setting bounds upon £. This rather unexpected effect is due to the amplifying factor (m/7). 

Another possibility for observable effects due to £ is to look at the averaged degree of circular polarization 

_ (p + h-p_h) 
110 - <p + M + P-H)' (43) 

where P\(w) is the total power distribution per unit frequency and polarization A, so that P\(lo) — P rn \/ujQ. The 
average here is calculated with respect to an energy distribution of the relativistic electrons, which we take to be 
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N{E)dE — CE p dE, in some energy range E\ < E < Ei. A standard choice for the exponent is 2 < p < 3. The 
result is 



n, 



fiU) 



n(p), 



(44) 



where Tl(p) will be estimated in what follows. Again, we have the presence of an amplifying factor in Eq. I|44f) . given 
by (fioj/qB), which is independent of the form of H(p) and not necessarily a small number. An estimation of this factor 
in the zeroth-order approximation (£ = = ff), which is appropriate in l|44|l . yields (ficu/qB) = lo/(ujqj) =771/7. The 
expression (|44|) is analogous to the well-known avera ge o f the degree of linear polarization Ulin = (p + l)/(p + 7/3), 
under the same energy distribution for the electrons |22j . 

Now let us compute the degree of circular polarization i|43|) in order to verify our previous statement. To this 
end we start from Eq. I|40[) and set 0=1 = n(Xz m ) everywhere, except in the critical terms involving rh\ c from 
where the corrections arise. We should also take into account that most of the radiation arises from the terms with 
m w rh\ c >> 1, where K 2 /3(l) — 0.49, «;(!) = 0.65. Thus the dominant term is 



Px(to) = D 



( 1 



\m\ c 



2/3 - 



K,{m/mxc) 



(45) 



where D is a constant. Starting from the definition (|43() . making an expansion in u = A£tt> and recalling that 
iti\c = rh\ c {u), we find 



n = fr 



du 



£=0, t/=0 



+ 0(C 2 , ^,...). 



(46) 



Since n(u) — y/l + u 2 + u we can rewrite d/du in terms of d/dn. Besides, fh\ c depends only on the combination nf3 
so that we can further go to d/d/3 obtaining 



dPx(to) 



du 



D 



' 2n 2 (3' 


d 


1 + n 2 n 





2/3 



-K(m/m c ) 



(47) 



where the first parenthesis on the RHS gives 1 upon evaluation, and we have already taken the limit n = l,fj = in 
the second parenthesis. In this way we now have 



T7 



E = M7. 



(48) 



Next, we consider the variable 



m 

TTLr 



2 LU 

3 UJQ 



1 -* = A 2 1 



A 2 = 



2 / fiuj 

3 \q~B 



(49) 



where A 2 is proportional to the amplifying factor in (|44|l . The above leads to 7 = Ax x l 2 . Again, a successive change 
of independent variables yields 



d/d/3 = /3 7 3 d/d-y = -2 A 2 d/da 



(50) 



where we have set f3 = 1 in the corresponding factor. Finally, we have to analyze the term to^ 3 ^ that multiplies 

k(x) inside the square bracket of Eqs. (|45|l and 147|) . Substituting m c = 37 3 /2 = 3 A 3 x~ 3 / 2 /2 we find that ^xm^ 3 ^ = 
(3A 3 /2) 2 / 3 is a constant that cancels when taking the ratio in l|46|) . Thus we are left with 



dK(x) 
dx 



(k(x)) ' 



(51) 



where the average (.) over E 1 has been replaced by one over x, via the change of variables E = 11 Ax 1 / 2 , with all 
the constant factors cancelling in the ratio of the two integrals. It should be pointed out that the factor A 2 arising 
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from the derivative in the numerator and responsible for the amplifying factor survives after taking this ratio. For 
the purpose of making an estimation of II (p) we take the energy range to be < E < oo. We then obtain 

4~ //xuA j°° x ^- 3 ^ 2 dK J x) dx 
110 = "3^ [WJ fcx( P -^Z)dx- (52) 



x» k(x) dx = ^-T + l)r + ~) , ,, + I >S - I . (.",:}) 



Using the expresionp^ 



o fx + 2 \2 3J \2 3 

and comparing with Eq. 1441) . we finally get 

(p-3)(3 P -i) (p + i)r(| + g)r, 
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n ^ = V %^T 7 r ; T^TTT W^mw^ll ■ P > 7 /3 (54) 



3(3p-7) (p-l)r(| 



12/ V 4 1 12/ 



The constraint p > 7/3 is required to avoid the divergence of the integral in the numerator of Eq. (|52|l at x — 0. On 
more realistic grounds, one should avoid the infinite upper limit of E, (the zero lower limit of x). In this case the 
divergence of the integral at x — 0, and hence the mathematical constraint on p, disappears, but the expression for 
n(p) becomes more complicated. In any case, the most important feature of the result in Eqs. (|44|l and H54|) is the 
presence of the amplifying factor (jjio/qB), which is independent of such details. 



IV. SUMMARY AND OUTLOOK 



We have presented a summary of a complete calculation of synchrotron radiation in the Myers-Pospelov effective 
model. The results can be understood in terms of the standard electrodynamics description of a parity-violating 
birrefringent media with a given helicity-dependent refraction index n(\£u) containing corrections arising from the 
electromagnetic sector of the theory. For each polarization-dependent refraction index, the results of Ref. |l7J are 
recovered. In particular we obtain that the opening angle of each contribution to the polarized radiation is given by 

59 S3 m- 1/3 « m~ c 1/3 = [1 - p 2 (E)n 2 (\iuj)} 1/2 . (55) 

Corrections also arise from the charged particle sector and are encoded in the modified dependence of the /3 factor 
upon the particle energy, according to Eq. (jSJ). 

In the case of the CRAB nebula, Tabled tells us that (,uj = ^oj/Mp « (r'/r) 2 , even for the upper bound of 
£ w 10~ 4 given by Gleiser and Kozameh in Ref. [§(]• In this way, consistency requires that the phase of the Green 
function l|39|) reduces to w (r — n ■!•'). This is equivalent to set £ = 0, n — 1 in all the arguments of the Bessel 
functions. This leads to 

Se^j-^E), Wc = ^ 7 3 (£), (56) 

where j(E) still includes corrections depending upon the fermionic parameter f) according to Eq. In others words, 
the corresponding results in Ref. are recovered for this particular situation. 

It would be interesting to search for astrophysical objects in which the amplifying factors are important, focusing 
on the observation of the average circular polarization, of possible corrections to the average linear polarization and 
of the Stokes parameters, among other possibilities. 
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